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Abstract 

We compute the static-light meson spectrum with Nf = 2 flavours of sea quarks using Wilson 
twisted mass lattice QCD. We consider five different values for the light quark mass correspond- 
ing to 300 MeV < mps < 600 MeV and we present results for angular momentum j = 1/2, j = 3/2 
and j = 5/2 and for parity V = + and V = —. We extrapolate our results to physical quark 
masses and make predictions regarding the spectrum of B and Bg mesons. 



1 Introduction 



A systematic way to study B and Bg mesons from first principles is with lattice QCD. Since 
arrih > 1 at currently available lattice spacings for large volume simulations, one needs to use 
for the b quark a formalism such as Heavy Quark Effective Theory (HQET) or Non-Relativistic 
QCD. Here we follow the HQET route, which enables all sources of systematic error to be 
controlled. 

In the static limit a heavy-light meson will be the "hydrogen atom" of QCD. Since in this limit 
there are no interactions involving the heavy quark spin, states are doubly degenerate, i.e. there 
is no hyperfine splitting. Therefore, it is common to label static-light mesons by parity V and 
total angular momentum of the light degrees of freedom j with j = |/ it 1/2|, where / denotes 
angular momentum and ±1/2 the spin of the light quark. An equivalent notation is given by 
S = (1/2)-, P_ = (1/2)+, P+ = (3/2)+, D_ = (3/2)-, ... The total angular momentum of 
the static-light meson is either J = j + 1/2 or J = j — 1/2, where both states are of the same 
mass. Note that in contrast to parity, charge conjugation is not a good quantum number, since 
static-light mesons are made from non-identical quarks. 

The static-light meson spectrum has been studied comprehensively by lattice methods in the 
quenched approximation with a rather coarse lattice spacing [1]. Lattice studies with Nf = 2 
flavours of dynamical sea quarks have also explored this spectrum [SJ [31 [U [U El [7] . Here (cf. also 
[8]) we use Nf = 2 and are able to reach lighter dynamical quark masses, which are closer to the 
physical u/d quark mass, so enabling a more reliable extrapolation. Note that in this formalism, 
mass differences in the heavy- light spectrum are 0{a) improved so that the continuum limit is 
more readily accessible. 

In this paper, we concentrate on the unitary sector, where valence quarks and sea quarks are of 
the same mass. This is appropriate for static-light mesons with a light quark, which is u/d. We 
also estimate masses of static-light mesons with light s quarks, albeit with a sea of two degenerate 
s instead of u and d. Within the twisted mass formalism, it is feasible to use Nf = 2 + 1 + 1 
flavours of dynamical sea quarks, which will give a more appropriate focus on the static-strange 
meson spectrum with light sea quarks. 

In HQET the leading order is just the static limit. The next correction will be of order l/mq, 
where mg is the mass of the heavy quark. This correction is expected be relatively small for 
b quarks, but larger for c quarks. Lattice methods to evaluate these I/tuq contributions to 
the B meson hyperfine splittings have been established and tested in quenched studies [HI llOj . 
We intend to explore these contributions using lattice techniques subsequently. An alternative 
way to predict the spectrum for B and Bs mesons is to interpolate between D and Dg states, 
where the experimental spectrum is rather well known, and the static limit obtained by lattice 
QCD assuming a dependence as l/mq. Thus the splittings among B and Bg mesons should be 
approximately rric/rrn, ~ 1/3 of those among the corresponding D and Dg mesons. 

For excited Dg mesons, experiment has shown that some of the states have very narrow decay 
widths [11]. This comes about, since the hadronic transitions to DK and DgM (where M is 
a flavour singlet mesonic system, e.g. rj' , vrvr or /o) are not allowed energetically. The isospin 
violating decay to Dgir together with electromagnetic decay to Dg^ are then responsible for the 
narrow width observed. A similar situation may exist for Bg decays and we investigate this 
here using our lattice mass determinations of the excited states. This will enable us to predict 
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whether narrow excited Bs mesons should be found. 

As weh as exploring this issue of great interest to experiment, we determine the excited state 
spectrum of static-light mesons as fully as possible. This will help the construction of phe- 
nomenological models and will shed light on questions such as, whether there is an inversion 
of the level ordering with /+ lighter than L at larger I or for radial excitations as has been 
predicted [El m H [E] . 

Since we measure the spectrum for a range of values of the bare quark mass parameter for the 
light quark, we could also compare with chiral effective Lagrangians appropriate to HQET. This 
comparison would be most appropriate applied to heavy-light decay constants in the continuum 
limit, so we will defer that discussion to a subsequent publication. 

This paper is organised as follows. In section [2] we review some basic properties of twisted mass 
lattice QCD. Moreover, we discuss particularities arising in static-light computations as well as 
automatic 0{a) improvement. In section [3] we present technical details regarding static-light 
meson creation operators and the corresponding correlation matrices we are using. We also 
explain how we extract the static-light spectrum from these correlation matrices and present 
numerical results for a range of light quark masses. We extrapolate these results both to the 
physical u/d quark mass and to the physical s quark mass. In section U we make predictions 
regarding the spectrum of B and Bg mesons by interpolating in the heavy quark mass to the 
physical b quark mass using experimental results as input. We close with a summary and a brief 
outlook (section [5]) . 



2 Twisted mass lattice QCD 
2.1 Simulation details 

We use X T = 24^ x 48 gauge configurations produced by the European Twisted Mass 
Collaboration (ETMC). The gauge action is the tree-level Symanzik (tlSym) action [16] 

Sg[U] = ^(boY.Tr(l-P'-\x-,f,,i^))+hY^TT[l-P'>^\x-,,,,u))\ (1) 

with the normalisation condition 5o = 1 — 85i and bi = —1/12. The fermionic action is the 
Wilson twisted mass (Wtm) action [I7l [181 [H] with Nf = 2 degenerate flavours 

Sf[x,X,U] = a'^^x{x)(^Dw + i^q'y5T3^x{x), (2) 

X 

where 

Dw = I (tm ( + - a V* V^) + mo , (3) 

V^j and V* are the standard gauge covariant forward and backward derivatives, mo and are 
the bare untwisted and twisted quark masses respectively and x = (x^"^ ) X^^^) represents the 
fermionic field in the so-called twisted basis. It is useful to introduce at this point the twist angle 
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io given by tan a; = /Ur/ttir, where //r, and rriR denote the renormahsed twisted and untwisted 
quark masses. This angle characterises the particular lattice action and must be kept fixed up 
to 0{a), while performing the continuum limit. 

The results presented in this paper have been obtained with gauge configurations computed at 
/5 = 3.9 corresponding to a lattice spacing a = 0.0855(5) fm. We consider five different values 
of fXq with rriQ tuned to its critical value at /iq = 0.0040 [201 \2T\ [22] (cf. Table [H where for each 
value the corresponding "pion mass" mps and number of gauge configurations is listed). With 
this tuning our target continuum theory is given by 

^ = xi^){lfiD^ + ifiK'y5r3^xix), (4) 

which is parameterised by the renormahsed twisted quark mass //r. The tuning guarantees 
automatic 0{a) improvement for physical correlation functions involving only light fermions 
|18| . In section [231 we will argue that automatic 0{a) improvement also holds for static-light 
spectral quantities without additional complications. 





mps in MeV 


number of gauge configurations 


0.0040 


314(2) 


1400 


0.0064 


391(1) 


1450 


0.0085 


448(1) 


1350 


0.0100 


485(1) 


900 


0.0150 


597(2) 


1000 



Table 1: bare twisted quark masses /iq, pion masses mps and number of gauge configurations. 



2.2 Static-light correlation functions 

To compute correctly a static- light correlation function with the Wtm lattice action ([2|), we 
follow the general procedure described in and reviewed in [19]. The procedure reads: 

(1) start with the continuum static-light correlation function you are interested in, 

(2) perform the axial rotation 

tp = ex.p (iu;-f5T3/2^x , V' = X^W (i^lbTs/^^ (5) 

on the fields appearing in the correlation function with a given value for lo, 

(3) compute the resulting correlation function with the Wtm lattice action with a choice 
of quark masses, such that tanw = /UR/mR up to 0{a), 

(4) perform the continuum limit with renormalisation constants computed in a massless scheme, 
tuning the untwisted bare quark mass in order to achieve the desired target continuum 
theory, i.e. the desired value of the twist angle lo. 
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Each value of uj defines a different discretisation, but when the continuum hmit is performed 
the result will be exactly the initially chosen static-light correlation function in the continuum 
with quark mass M"^ = m^^ + /x^. 

In the following we give an explicit example. In QCD the pseudoscalar and scalar static-light 
currents read 

V^^\x) = Q{x)-f5iJ^''\x) , S'^^'^x) = (6) 

where Q is the static quark fielc|3 and tp^^^ is a single flavour of the light fermion doublet 
iP = (tpi^) , ^W). Let us suppose we are interested in computing in continuum QCD the static- 
light pseudoscalar-pseudoscalar correlation function 

Cvv = ((^^*^*)R(x)(P^*^^)J,,(y))^^^^_^^, (7) 

where we write an index (Mr,0) to specify that the continuum action has a vanishing twisted 
mass and a renormalised untwisted mass given by Mr. We perform the axial rotation ([5]) 
obtaining 

cos2(w/2)Z^Cpp + sin2(u;/2)Z|Css - i cos(w/2) sin(tj/2)ZpZs (Cps - Csp) , (8) 

where Zp and Zs are the standard renormalisation constants for static-light currents computed 
in a massless scheme with Wilson fermions. Note that for the static-light case, Zy = Zp and 
Za = Zs- This correlation function has to be computed with the Wtm action ([2]) with quark 
masses tuned accordingly to the value of oj chosen. The Cxx correlation functions in ^ are 
defined in terms of currents in the twisted basis 

Cpp = (p^'-\x){P^'-')\y)) Css = (s^'^\x){S^'^')Hy)) ^ , . . . , (9) 

where 

= Q-/5x'-"'\x) , 5"*^* (a;) = Q{x)x^''\x). (10) 

Once the continuum limit of the correlation function ([8]) has been performed, the result will be 
the original correlation function ([7]) with M-^ = rn^ + fi^. 

However, to compute spectral quantities it is sufficient to analyze a matrix of correlation func- 
tions of bare currents with the appropriate quantum numbers. We will discuss this in detail in 
section [22 



2.3 Automatic 0{a) improvement of static-light meson masses 

Spectral quantities like hadron masses extracted from lattice simulations of Wilson fermions will 
in general be affected by 0{a) discretisation errors. In the particular case of masses extracted 



^We will discuss the static quark action in section [3. 2. II 
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from static-light correlation functions the 0{a) discretisation errors come from the dimension- 
5-operators of the Symanzik effective action of the light and static quarks. 

The Symanzik effective action for the Eichten-Hill (EH) static action contains only one term, 
which contributes to the 0{a) corrections of the linearly divergent static self-energy |23j . In this 
paper all observables we consider are differences, where this static self-energy cancels. Moreover, 
this result is independent on the particular lattice static action chosen, as long as it preserves 
the relevant symmetries of the EH action. This is the case for our choice of static action (cf. 
section 13.2. 1|) . 

As a consequence, the only 0{a) errors which could affect our results, come from the dimension- 
5-operators of the Symanzik effective action of the light quarks. The light quark action used in 
this paper is Wtm at maximal twist. It is by now well known that at maximal twist a single 
insertion of a dimension-5-operator of the Symanzik effective action into parity even correlation 
functions vanishes, because, independently on the lattice basis adopted, these operators are 
parity odd and the insertions have to be evaluated in the continuum theory, where parity is 
a preserved symmetry [H]. We can conclude that all the spectral quantities, when the static 
self-energy has been removed, are automatically 0{a) improved. 



2.4 Spectral decomposition and parity mixing 

In this section we explain, how to analyze lattice results for static-light correlation functions 
obtained in the twisted basis. In particular we concentrate on the assignment of parity labels 
to extracted static-light meson states. 

We start from the physical basis and, for simplicity, consider only two operators, the pseudoscalar 
and the scalar static-light current, and only two states, which we label by |1) and |2). The 
explanation carries over to the more general case in a straightforward way. 

First consider the following matrix of correlation functions in the physical basis: 

- ( S(') ) ■ 

where C'p'p{t) has been defined in ([7]) with x = (t,0) and y = (0,0) and analogously the others. 
The parity of the operators (7^'^*'^'')r and (5'^*'^'')r is determined by the parity transformation 
properties of the associated field, i.e. {V^^^^)k has negative parity and (5'^*'^')r has positive 
parity. Even if parity is broken at finite lattice spacing, one can still assign a parity label to 
each of the states we use to decompose the correlation functions [18] . If we consider only two 
states, the spectral decomposition will have the form 

CiT^ - ( 

^ ' ~ V («f)*«r i«fp 

where we have defined 

{al^r = {^\V''-'\l,2) , {al,r = {n\S^'-'\l,2). (13) 

The correlation functions C-ps and Cgp vanish in the continuum limit, because parity is a sym- 
metry of QCD. This means by universality that at finite lattice spacing they are at most of 



1 ^-A/it , / I {0-2 )*a2 \ -M2t 
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0{a). Since C-p-p and C55 are of 0{1) in the continuum limit, we can conclude that for given n 
either is of 0{1) and af is of 0{a) or the opposite way round [18]. We can conclude that if 
is of 0(1), the state |n) has the same parity as the formal parity of 'P^^'^^^ which in this case 
is negative. Moreover, is of 0{a) and has to vanish in the continuum limit. 

We now perform the axial transformation ([5]). The relation between correlation functions up to 
discretisation errors is, for example, for C-p-p 



Cpp = cos^(a;/2)Z^Cpp + sin^(u;/2)Z|Css -«cos(w/2)sin(w/2)ZpZs(^Cps -Cspj. (14) 
For the matrix of correlation functions in the twisted basis 

- ( csp'S ) <^^> 

we can also perform a spectral decomposition considering again only the states |1) and |2): 



From ()12p . (jl4p and (jl6p we can conclude 



\a 



= cos2(cj/2)Z^|6f_2l^ + sin^('^/2)^s|fei,2p + 2cos(w/2)sin(tj/2)ZpZsIm((6f_2)*fci,2j 

(17) 

af2p = cos2(u;/2)Z||6f 2!^ +sin^(^^/2)^p|^r2p + 2cos(cj/2)sin(cj/2)ZpZsIm((6f 2)*fci'2^ 



(18) 

If the state |1) has negative parity, |af p has to vanish as 0{a?') in the continuum limit, while 
|a^P has to be of 0{\). Since the first two terms on the right hand side of (jl7p are positive and 
non-vanishing in the continuum limit, there must be a cancellation coming from the third term. 
In fact we immediately see that this third term has opposite sign for |a^2p compared to |af 2p- 
This allows us to identify the parity of the states |1) and |2) without knowing the exact values 
of the renormalisation constants and the twist angle. The criterion will be the following: if 



Im 



{(hlTht) < 0, (19) 



the state |1) has negative parity, otherwise positive parity. The other cases follow accordingly. 

This method, which we have described for a simple case, is valid independently of the number 
of states considered and the kind of operators studied. At finite lattice spacing it provides a 
way to assign a formal parity to each of the extracted states. 

The method extends to all cases, where the light degrees of freedom involve fermions in the 
twisted basis, e.g. for static-light mesons, but also for baryons. 
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3 The static-light meson spectrum 



3.1 Static-light trial states 

3.1.1 Static-light meson creation operators in the continuum 

It is convenient to discuss static-light mesons treating the static quark as a four component 
spinor since the symmetries of hadronic bilinears are well studied [23]. In the continuum an 
operator creating a static-light meson with well defined quantum numbers J , j and V is given 
by 

0^^\x.) = (5(x) j dhr{ii)U{yL;x + rn)ip^''\x + rn). (20) 

(5(x) represents an infinitely heavy antiquark (here a Dirac spinor) at position x, J dh denotes 
an integration over the unit sphere, C/ is a straight parallel transporter and ip^^^x + rh) creates 
a light quark at position x -|- rh separated by a distance r from the antiquark (of course, using 
ijjid-) instead of V'^"'' would yield identical results). F is an appropriate combination of spherical 
harmonics and 7 matrices coupling angular momentum and quark spin to yield well defined 
total angular momentum J (static quark spin included) and j (static quark spin not included) 
and parity V. The meson creation operators used in the following are listed in Table [2j 



r(h) 




f 


Oh 


lattice 


notation 






(1/2)- 


Ai 


(1/2)- , (7/2)- , ... 


S 


1 , -fjhj 


0+ [1+] 


(1/2)+ 




(1/2)+ , (7/2)+ , ... 


P_ 


71^1 — 72^2 (and cyclic) 


2+ [1+] 


(3/2)+ 


E 


(3/2)+ , (5/2)+ , ... 


P+ 


75(71^^1 - 72^-2) (and cychc) 


2- 


(3/2)- 




(3/2)- , (5/2)- , ... 


D± 


7l^2T^3 + 72%"! + 73nin2 




(5/2)- 




(5/2)- , (7/2)- , ... 


D+ 


75 (71 "-2% + 72n3ni + 73nin2) 


3+ [2+] 


(5/2)+ 




(5/2)+ , (7/2)+ , ... 


F± 



Table 2: Static-light meson creation operators. The other mesonic states that are degenerate 
with that created are noted in square brackets. 



3.1.2 Static-light meson creation operators on a lattice 

Here we present the construction of appropriate lattice operators to create the states of interest, 
following [H [231 . When putting static-light meson creation operators (|20p on a lattice, one 
has to replace the integration over the unit sphere by a discrete sum over lattice sites, which 
have the same distance from the static antiquark at position x. For the operators in Ai and E 
representations we use six lattice sites, i.e. 

0(^)(x) = Q(x) F(h)C/(x;x + rn)x^"^(x + rn), (21) 

n=±ei,±e2,±e3 
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whereas for those in the A2 representation one has to use eight lattice sites, i.e. 

0(r)(x) = Q(x) r{n)U{x;x + rn)x^''\x + rn). (22) 

n=±ei±e2±e3 

In the first case the spatial parallel transporters are straight paths of links, while in the second 
case we use "diagonal links", which are averages over the six possible paths around a cube 
between opposite corners projected back to SU(3). 

The states created by these lattice meson creation operators do not form irreducible represen- 
tations of the rotation group SO (3), but of its cubic subgroup Oh. Therefore, these states have 
no well defined total angular momentum, but are linear superpositions of an infinite number 
of total angular momentum eigenstates. The common notation of the corresponding Oh rep- 
resentations together with their lowest angular momentum content are also listed in Table [21 
Note that we do not consider Oh representations Ti and T2, because these representations 
yield correlation functions, which are numerically identical to those listed (e.g. Ti would be 
r = 7j or r = 757j, which gives the same correlations as F = 75 and T = 1, and T2 would be 
r = 7in2 + 72?^! or r = 75(71712 + 72711), which gives the same correlations as T = 7irei — 72^12 
and r = 75(7ini - 72^2)). 

Since the D_ and the states as well as the and states are expected to have a 
similar mass, we do not have unambiguous lattice operators to determine D_ and F_ but rather 
operators, which have an admixture of D± and F± respectively. We label these operators as D± 
and F± (cf. Table [2]). 

We have also replaced the light quark fields in the physical basis V''-"^ by their counterparts in 
the twisted basis x^"''- Note that trial states created by such twisted basis operators are not 
eigenstates of parity. Nevertheless, as we have discussed in section 12. 4^ it is possible to assign 
unambiguously a parity label to the masses extracted from the time dependence of such twisted 
basis correlators. 

3.1.3 Smearing techniques 

When performing a lattice study of the static-light meson spectrum, the following points have 
to be considered: 

• It is imperative to use trial states with large overlap to low lying energy eigenstates. Only 
then the corresponding meson masses can be extracted from correlation functions at small 
temporal separations, where signal-to-noise ratios are acceptable. 

• To determine excited states for a given Oh representation, it is necessary to have a whole 
set of linearly independent trial states belonging to that Oh representation. 

To fulfill both requirements we use different "radii" r (cf. eqns. (12ip and (122p ) and apply APE 
smearing and Gaussian smearing also with different parameters. The resulting extended trial 
states have significantly better overlap to low lying energy eigenstates than their unsmeared 
counterparts. 
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APE smearing of spatial links 

After Aape iterations APE smeared spatial links [25] are given by 

^ i=±l,±2,±3 



(23) 



where U^^^ are the original unsmeared links, gape is a weight parameter and Psu(3) denotes a 
projection back to SU(3) defined by 

with det(i7')^/'^ being that root closest to 1. 
Gaussian smearing of light quark operators 

After Acauss iterations Gaussian smeared light quark operators [261 [27] are given by 

^(-'VGauss)^^;) = 



1 



1 + 6k 



j=dzl,±2,±3 ^ 



where x^^^ are the original unsmeared light quark operators and [/(^ape) denote APE smeared 
spatial links. 

3.2 Correlation matrices 

For each Oh representation we compute 6x6 correlation matrices 

CKK'it) = /O(^)(t)(0(^'))t(0)\, (26) 



where O(^) is a static-light meson creation operator (cf. eqns. (j2ip and (j22p ) with K denoting its 
parameters, i.e. K = {T , A^auss j f) (we have chosen Aape = 10, oape = 0.5 and ^causs = 0.5 
for all operators). Detailed information about the operator content of the correlation matrices 
is given in Table [3l 



The width of a Gaussian smeared light quark operator (j25p in lattice units is approximately 
given by 



/2Ag auss'^'Gauss 
V 1 + 6KGauss 
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Oh 


r 


-^Gauss 


r 




R in fm 


A-. 


75 


ou 


Q 


^ fil 


n 4s 






DU 


D 


Q nn 
y.uu 


n 77 




1 


30 


3 


5.61 


0.48 






fin 

uu 


fi 


Q nn 


n 77 






30 


3 


5.61 


0.48 




7i^j 


30 


3 


5.61 


0.48 


E 


71X1 — 72^2 (and cyclic) 


30 


3 


5.61 


0.48 






60 


6 


9.00 


0.77 






90 


3 


8.74 


0.75 




75(7ia;i — 72a;2) (and cyclic) 


30 


3 


5.61 


0.48 






60 


6 


9.00 


0.77 






90 


3 


8.74 


0.75 


A2 


+ 72X3X1 + 73X1X2 


30 


2 


5.88 


0.50 






60 


4 


9.64 


0.82 






90 


2 


8.91 


0.76 




75(71 a;22;3 + 72^^32:1 + 73^:1x2) 


30 


2 


5.88 


0.50 






60 


4 


9.64 


0.82 






90 


2 


8.91 


0.76 



Table 3: static-light meson creation operators used for the Ai, E and A2 correlation matrices. 



For KGauss = 0.5 and A^causs = (30 , 60 , 90) this amounts to a ^ (2.74 , 3.87 , 4.74). Taking also 
the parameter r into account one can estimate the radius of a static-light trial state: 
R/a = Vr"^ + 3(7^ for the Ai and E representations and R/a = VSr"^ + 3cr^ for the A2 represen- 
tation. The radii of the trial states used are also listed in Table [3] both in lattice units and in 
physical units. 

Note that to identify the parity of states extracted via fitting it is important to compute corre- 
lation matrices, which contain for each operator T also its counterpart 75r (cf. section [27i|) . 

3.2.1 Quark propagators 

When evaluating the correlations (p6]) . both static quark propagators and light quark propaga- 
tors appear. To improve signal-to-noise ratios, we apply the following techniques. 
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Static quark propagators 

To improve the signal to noise ratio for static-light correlation functions, we use the HYP2 static 
action [28l [291 EO] • Static quark propagators are given by 

{Q{x)Q{y))^__ = 5^^\-K-y)U^''^''''\x-y)(^Q{y,-xo)'^^ + e{xo-yo)^^y (28) 

where (. . .)q q denotes the integration over the static quark fields and U{x; y) is a path ordered 
product of HYP2 smeared links along the straight path from x to y. 



Light quark propagators 



To exploit translational invariance, it is imperative to use stochastic methods for the light quark 
propagators. The correlators can then be evaluated at a large number of source points, while 
only a few inversions of the lattice Dirac operator have to be performed. One very powerful 
method is maximal variance reduction [1]. A somewhat easier method to implement is to use 
stochastic sources on time slices and this has been found to give reasonable results |31j . Because 
we have inverted from such time-slice sources as part of our light-light meson studies [20^ [2H [22] , 
we follow this latter route, since it is computationally much quicker for us. 

For each gauge configuration we use Ng stochastic Z2 x Z2 sources a = l,...,iV^ located 
on the same timeslice. For our lightest three /ig values we take A^^ = 4 sources, which are the 
same for each of the four spin components so that we can re-use previous inversions [201 121[ [22]. 
For our heavier two values, we had to redo the inversions so we use only Ns = 1 source with 
random values in each of the spin components. 

After solving 

D^^\^{x-y)<l^^^\y) = C^^\x), (29) 

where D^l^ = Dy^ + ifi^fi^ is the twisted mass Dirac operator acting on x^"\ the light quark 
propagator is given by the unbiased estimate 



(x^"^(x)x(")(y)) _ = {D^^\j-\x-y) « j;<A(")(x)(e("))t(y), (30) 
where (. . ■)x,x denotes the integration over the light quark fields. 



3.3 Extracting static-light meson masses from correlation matrices 

Assuming that for sufficiently large t the correlation matrix ([26]) can be approximated by the n 
lowest lying energy eigenstates i = 1, . . . , n we use the ansatz 



K\ 



I ■ 



(31) 



1=1 
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The correlation matrix (j26p in terms of the ansatz is 



CKK'{t) 



Eit 



(32) 



The parameters Ei and bf are determined by minimising 

2 _ [ CKK'jt) - CKK'{t) \ /„oN 

- ^ ^ 1 a{CKK'{t)) 1 ' ^^^^ 



where (j{CKK'{t)) denotes the statistical error of CxK'it)- 

In the following we apply this fitting procedure with n = 4 exponentials. To obtain physically 
meaningful results with small statistical errors, it is essential to determine an appropriate fitting 
range tmin ■ • • ^max- To this end, we have performed correlated fits with various fitting ranges 
using eigenvalue smoothed covariance matrices |32j . We have found that tmin = 3 gives reason- 
able reduced values (cf. Tabled]), while data points beyond tmax = 12 seem to be dominated 
by statistical noise, i.e. including them in the fits does not alter resulting meson masses nor 
corresponding statistical errors. 



Oh 


I2q = 0.0040 


/iq = 0.0064 


= 0.0085 


flq = 0.0100 


flq = 0.0150 


Ai 


1.89 


2.30 


2.35 


0.95 


1.16 


E 


1.21 


1.33 


1.70 


2.04 


2.09 


A2 


1.56 


1.96 


1.28 


1.16 


1.26 



Table 4: x^/dof from correlated fits for different Oh representations and different /iq. 



As has already been discussed in section 13.1.21 it is difficult to unambiguously determine the 
total angular momentum j of a state obtained from a lattice computation. This is, because 
for every Oh representation there exists an infinite number of possible total angular momentum 
eigenstates (cf . Table [2]) . In the following, we assume that the low lying states we are going to 
study have the lowest total angular momentum possible, i.e. we assign j = 1/2 to states from 
Ai, j = 3/2 to states from E and j = 5/2 to states from A2. Parity on the other hand can 
directly be read off from the coefficients bf (cf. section [2^ . 

Since static-light meson masses diverge in the continuum limit due to the self energy of the 
static quark, we always consider mass differences, where this self energy cancels. Mass differences 
between various static-light mesons with quantum numbers and the lightest static-light meson 
((l/2)~ = S ground state) for all five /iq values are collected in Figure [Hand Table [5l Statistical 
errors have been computed from 100 bootstrap samples. 

To check the stability of the fitting method, we have performed computations with different 
parameters (number of states n, fitting range ^min • • • ^max; operator content of the correlation 
matrices). We have obtained results which are consistent within statistical errors. 
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3 


/iq — U.UU4U 


flq — U.UU04 


/iq — U.UUoO 


/iq — U.UIUU 


/iq — U.UloU 


(1/2)"'* = S* 
(1/2)+ ^ P_ 


777(17) 
389(16) 


808(19) 
428(12) 


839(22) 
447(10) 


780(34) 
456(17) 


782(32) 
495(16) 


(3/2)+ ^ P+ 
(3/2)- = D± 


473(10) 
813(24) 


496(8) 
828(19) 


488(7) 
833(16) 


486(12) 
861(27) 


479(14) 
858(21) 


(5/2)- ^ 
(5/2)+ ^ F± 


823(24) 
1134(35) 


887(14) 
1205(27) 


887(15) 
1173(24) 


862(24) 
1136(34) 


846(42) 
1205(28) 



Table 5: static-light mass differences m(j^) — 7n[S) in MeV for different ^q. 
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Figure 1: static-light mass differences linearly extrapolated to the physical u/d quark mass and 
the physical s quark mass. 



3.4 Extrapolation to physical light quark masses 

We linearly extrapolate our static-light mass differences in (mps)^ to the physical u/d quark 
mass (mps = 135 MeV) and the physical s quark mass (taken here as mpg = 700 MeV). Results 
are shown in Figured] and Table [6l We also list the corresponding x^/dof values indicating that 
straight lines are acceptable for extrapolation. A more thorough study using extrapolations 
based on chiral effective theories will be attempted when we are able to extract the continuum 
limit of our results at each light quark mass value. 

Note that we consider the unitary sector, where valence quarks and sea quarks are of the same 
mass. For the s quark extrapolated results this implies a sea of two degenerate s instead of a 
sea of u and d. If the sea-quark mass dependence of our spectra is small, as usually assumed, 
then our results will be a good estimate of the physical static-strange meson spectrum. This 
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u/d quark extrapolation: 


s quark extrapolation: 




J 




?n(j^) — 771(5*) in MeV 


TTii^j^^ — mi^S") in MeV 


X /aoi 


(1/2)"'* ^ 




791(23) 


816(43) 


1.82 


(1/2)+ ^ 


P_ 


371(16) 


554(23) 


0.44 


(3/2)+ ^ 


^+ 


487(11) 


486(19) 


1.22 


(3/2)- ^ 




804(23) 


887(33) 


0.21 


(5/2)" ^ 


^+ 


864(27) 


894(50) 


2.24 


(5/2)+ ^ 




1149(33) 


1215(44) 


1.40 



Table 6: static-light mass differences linearly extrapolated to the physical u/d quark mass and 
the physical s quark mass. 



limitation can be removed, in principle, by performing similar computations on A'^^ = 2 + 1 + 1 
flavour gauge configurations, which are currently being produced by ETMC [33j . 

We have performed a similar extrapolation for the mass difference of the P wave states. When 
extrapolating to the physical u/d quark mass, we find m(P+) — m{P-) = 117(17) MeV, i.e. the 
P_ = (1/2)+ state is lighter than the P-f = (3/2)+ as usually expected. When increasing the 
mass of the light quark, we observe a reversal of this level ordering, m(P_)— m(P+) = 71(23) MeV 
at the physical s quark mass. It will be interesting to study this in the continuum limit, in 
particular since such a reversal is predicted by certain phenomenological models jl2 t ll3 1 [T ^ ll5j. 

In principle, our excited states could be two-particle states since we have dynamical sea quarks. 
In practice, the two-particle state is expected to be weakly coupled to the operators we use 
(which are constructed assuming one particle states). Some exploration of transitions to two 
particle static-light mesons has been made which confirms this expectation |31j . 

4 Predictions for B and Bg mesons 

To make predictions regarding the spectrum of B and Bg mesons, we interpolate between the 
static-light lattice results obtained in the previous section and experimental results for charmed 
mesonqj [11]. To this end, we assume a linear dependence in l/mq, where mq is the mass 
of the heavy quark. This interpolation introduces a possible systematic error, which, however, 
we consider to be smaller than the systematic errors coming from the continuum limit, the 
extrapolation to light quarks and the treatment of the strange sea. The most important of these 
systematic errors is that involved in the continuum limit and that will be quantified when we 
have results at finer lattice spacings. 

■^For the states B, D, D* , Dq and DJ experimental results for charged as well as for uncharged mesons exist. 
We use the average in the following. 
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4.1 B mesons 

Results of the interpolation between our u/d extrapolated P wave lattice results and experimen- 
tal results on D mesons are shown in Figure [2^ and Table El 

• To predict m{BQ) — m{B) and m{B\) — m{B), we interpolate between our static spin 
degenerate P- = (1/2)+ state, i.e. m{P-) — m{S), and experimental data on 

m{D^) - m{D) and m(Di(2430)°) - m{D). 

• To predict m{Bi) — m{B) and m(i?2) — m{B), we interpolate between our static spin 
degenerate P+ = (3/2)"'" state, i.e. m{P^) — m{S), and experimental data on 
m{Di{2A20)^) — m{D) and 772(1)2) ~ 'm{D). Here we assign the D\ states assuming that 
states with similar widths belong to the same multiplet. 

• The line labeled "S" = (l/2)~" in Figure[2^ shows that m{B*) —m{B) is lighter by a factor 
of ~ rric/rrH) than m[D*) — m{D) indicating that a straight line is a suitable ansatz for 
interpolation and that the estimate of mdm^, = 0.3 [TTj is reasonable. 

• A comparison with experimental results from CDF and D0 [34^135) on m{Bi) — m{B) and 
777(^2) ~ m[B) shows that our lattice results are larger by ~ 10% (cf. Table [7]). There 
is another resonance listed in [11] with unknown quantum numbers J'^, m{B*j) — 'm{B), 
which is rather close to our m{BQ) — m{B) and m{B^) — m{B) results. For a conclusive 
comparison it will be necessary to study the continuum limit, which will be part of an 
upcoming publication. 





m — m[B) in MeV 




777 — m{Bs) in MeV 


state 


lattice 


CDF 


D0 


PDG 


state 


lattice 


CDF 


D0 


PDG 


^0 


413(19) 








BsO 


493(16) 








Bl 


428(19) 








B*si 


535(16) 








Bi 


508(8) 


454(5) 


441(4) 




Bsl 


510(13) 


463(1) 








519(8) 


458(6) 


467(4) 




Bs2 


521(13) 


473(1) 


473(2) 




B} 








418(8) 


BIj 








487(16) 



Table 7: lattice and experimental results for P wave B and Bg states. Errors on lattice results 
are statistical only. 



4.2 Bg mesons 

For Bg mesons we proceed in the same way as for B mesons, using our s quark extrapolated 
static-light lattice results and experimental results on Dg mesons (cf. Figure [2}3 and Table [7]). 

• To predict m{B*Q) — m{Bg) and m{B*T^ — m{Bg), we interpolate between our static spin 
degenerate P_ = (1/2)+ state, i.e. m{P-) — m{S), and experimental data on 
m{D*o) - m{Dg) and m(L>^i(2460)) - m{Dg). 
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a) Excited B meson masses 
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b) Excited B meson masses 
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Figure 2: Static-light mass differences linearly interpolated to the physical h quark mass. 



• To predict m{Bsi) — m{Bs) and m{B*2) — m{Bs), we interpolate between our static spin 
degenerate = (3/2)+ state, i.e. m(P+) — m(S'), and experimental data on 
m{Dsi{2536)) — ni{Ds) and m{Ds2) — m{Ds). This time we assign the Dgi states according 
to the expectation that the splitting between Z)si("j = 3/2") and Ds2 is roughly 
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mb/nic ~ 3.3 times larger than that between Bgi and B*2, which is according to |36| 
approximately lOMeV. We also illustrate the opposite assignment in Figure [2)3 for com- 
pleteness. 

• The line labeled "5 = (1/2)"" in Figure [2)3 shows that m{B*) — m{Bs) is lighter by a 
factor of ~ rric/mf, than m{D*) —m{Ds) indicating that a straight line is a suitable ansatz 
for interpolation and that the estimate of mc/rrih = 0.3 |llj is reasonable. 

• A comparison with experimental results from CDF and D0 |36^ [37] on m{Bi) — m{B^) 
and m(i?|) — m{B^) shows that our lattice results are larger by ~ 10% (cf. Table[7|). There 
is another resonance listed in [11] with unknown quantum numbers J^, m{B*j) — ni{Bs), 
which is rather close to our m{B*Q) — m[Bs) result. For a conclusive comparison it will be 
necessary to study the continuum limit, which will be part of an upcoming publication. 

• We also plot the BK and B*K thresholds in Figure [2)3. The fact that our lattice results 
on the P wave states B*q, B*^, Bgi and B*2 are larger indicates that corresponding decays 
are energetically allowed. Therefore, one should expect that these states may have a larger 
width compared to the corresponding excited Dg states. 

5 Conclusions 

We have explored the low lying static-light meson spectrum using Nj = 2 flavours of sea quarks 
with Wtm lattice QCD. We have presented results for total angular momentum of the light 
degrees of freedom j = 1/2, j = 3/2 and j = 5/2 and for parity V = + and V = —. The lattice 
spacing is a = 0.0855(5) fm and we have considered five different values for the light quark mass 
corresponding to 300 MeV < mps < 600 MeV. 

We have extrapolated our results in (mps)^ both to the physical u/d quark mass and to the 
physical s quark mass. Moreover, we used experimental results from D and Dg mesons to 
interpolate in the heavy quark mass from the static case to the physical h quark mass. We are 
able to predict the spectrum of excited B and Bg mesons from first principles. Our formalism 
has lattice artifacts of order a? and we shall be able to control these in future work by studying 
smaller a values. Comparing our current predictions to available experimental results, we find 
agreement up to 10% with P wave B and Bg mesons. 

Throughout this paper we have considered the unitary sector, where valence quarks and sea 
quarks are of the same mass. Particularly for our Bg results, this implies a sea of two degenerate 
s instead of a sea of u and d. We plan to improve this by performing similar computations on 
Nf = 2 + 1 + 1 fiavour gauge configurations, which are currently produced by ETMC. Another 
important issue in the near future will be an investigation of the continuum limit, which amounts 
to considering other values for the lattice spacing. Such a study will be necessary for a conclusive 
comparison between lattice results and experimental results for B and Bg mesons. We also plan 
to compute static-light decay constants and to make a detailed comparison with chiral effective 
Lagrangians. 
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